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ERNST'S EQUATION FOR COLLIDING 
GRAVITATIONAL WAVES 

The solutions being considered here for colliding plane waves all have a 
pair of commuting Killing vectors that are assumed to exist globally. It 
may therefore be expected that the solutions obtained may be related to 
known cylindrically symmetric solutions, or to stationary axisymmetric 
solutions, which similarly have a pair of Killing vectors. Such a relation 
was first pointed out by Kinnersley (1975), and by Fisher (1980). The ex- 
act relation between these solutions has been established more recently by 
Chandrasekhar and Ferrari (1984), and Chandrasekhar and Xanthopoulos 
(1985 a), and exploited by these authors and their colleagues. 

In this chapter we will present an analysis of the colliding wave prob- 
lem using a method that has become familiar in the study of stationary 
axisymmetric space-times. In this case the field is described in terms of 
a complex potential function that is referred to as the Ernst potential 
(see Ernst 1968a). This approach leads directly to various methods for 
generating further exact solutions. 

11.1 A derivation of the Ernst equation 

First it may be recalled that, of the vacuum field equations (6.22a-/) 
considered in previous chapters, (6.22a) may immediately be integrated 
to give 

e~ u = /(«) + g(v), (11.1) 

and (6.22a 1 , e) are integrability conditions for the remaining equations. 
Attention is thus focused on these main equations for the metric functions 
V(u,v) and W(u,v). 

It is now appropriate to consider a different combination of these 
functions, by putting 

x(u,v) = e~ v sechW, u(u,v) = e~ v tanhW, (H-2) 

or, inversely 

e v = (u 2 + x 2 )" 1/2 , smhW = u/ X - (11.3) 
95 



96 Ernst's equation for colliding gravitational waves 

This modifies the form of the line element (6.20), which now becomes 

ds 2 = 2e~ M dudv - e~ u (x dy 2 + x _1 (dx - udy) 2 ) . (11.4) 
It is also convenient to introduce the complex function 

Z = X + iu. (11.5) 
With this, the line element (11.4) can be written in the form 

2e~ u 

ds 2 = 2e~ M dudv - — ^(dx + iZdy)(dx - iZ dy) (11.6) 

[Z + Z) 

and the two main equations (6.22d, e) can be written as the single complex 
equation 

(Z + Z)(2Z UV - U U Z V - U V Z U ) = AZ U Z V . (11.7) 

It may be seen that this is in fact Ernst's equation, which can be written 
in the coordinate-invariant form 

(Z + Z)V 2 Z = 2(VZ) 2 (11.8) 

where (V(j)) 2 = g^ v 4> ,v is the square of the gradient of an arbitrary 
scalar field 4> which, in this CclSG, IS ct function of the two (null) coordinates 
only. Similarly, V 2 is the 3+1-dimensional Laplacian operator (or the 
generalized d'Alembertian) given by 

Using the above notation, the non-zero components of the Weyl ten- 
sor given by (6.23) can be written as 

n = z J z +z)2 t {z + 2){Zvv ~ UvZv + MvZv) ~ 2Z ^ 

UuUy_ Z U Z 

2 4 (Z + Z) 2 K ' 

n = z J z +z)2 t {z + z){Zuu ~ UuZu + MuZu) ~ 2Z 3 • 

When considering Ernst's equation, it is frequently found to be useful 
also to introduce an associated function E, defined by 

1 +E Z - 1 , 

Z= -, or E=- . (11.11) 

1-E y Z+l v ; 
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With this, the line element (11.4) or (11.6) can be written in the alterna- 
tive form 

ds 2 = 2e~ M &u&v 

- yt^e- [t 1 - E ^ dx + ^ + E ^y\ K 1 - E ^ dx - ^ + E ^y\ 

(11.12) 

and the main equations (6.22d, e), or alternatively (11.7), become 

(1 - EE)(2E UV - U U E V - U V E U ) = —AEE U E V . (11.13) 

This is the alternative form of the Ernst equation and may be rewritten 
in the coordinate-invariant form 

(1 - EE)V 2 E = -2E(VE) 2 . (H-14) 

Using this function, the non-zero components of the Weyl tensor 
(11.10) can be written as 

U V E V +M V E V ) + 2EE 2 ) 

(11.15) 

U u E u + M u E u )+2EEl). 

The problem now involves finding appropriate solutions of (11.7) or 
(11.13). These equations, however, contain the derivatives of U, and hence 
they depend on the arbitrary functions f(u) and g(y) that are specified 
by the incoming waves. This apparently explicit dependence on initial 
conditions may in fact be removed by adapting the coordinate system of 
Chandrasekhar and Ferrari (1984). Here we again use (10.9-12) and put 

f(u) = i cos(V> + A), g(v) = \ cos(V> - A) (11.16) 

where ip and A are considered as time-like and space-like coordinates, 
which may then be rescaled by putting 



T — E 2 1 

»g = i/ i-g ( i- g g). (d-gg)(g~- 



11 11 EE 

to '~ J u'- J v J - J u J -'v 



4 (1-EE) 



1 - E 2 1 



t = simp, z = sin A. 



(11.17) 



98 



Ernst's equation for colliding gravitational waves 



In this coordinate system, equation (11.7) takes the more familiar explicit 
form of Ernst's equation: 



(Z + Z)[((l-t 2 )Z t ) t -((l-z 2 )Z z ) 



(11.18) 

= 2((l-t 2 )Z 2 -(l-z 2 )Z 2 ) 



and (11.13) similarly becomes 



= -2E ((1 - t 2 )E 2 t - (1 - z 2 )E 2 z ) . 

(11.19) 

The intermediate steps in the derivation of these equations may be de- 
duced from (16.6) and (16.11). See also (12.30-31). 

It may be noticed that, in this case, the solution of these equa- 
tions immediately determines some of the metric functions. This is in 
marked contrast to their application in stationary axisymmetric space- 
times, where the Ernst equation only determines potentials for the fields. 

The original metric functions, as considered in previous chapters, are 
now given by 

2V = = (i-g)(i-g) 

_< 1 + £ »< 1 + £ » (11.20) 

sinhVr = —i- =f = —%- 

(Z + Z) (1 - EE) 



11.2 Boundary conditions 

When looking for solutions of the Ernst equation for stationary axisym- 
metric space-times, it is appropriate to require that solutions be asymp- 
totically flat. However, for colliding plane waves very different boundary 
conditions apply. 

For colliding plane waves it is necessary to choose Z or E, and hence 
V and W, such that the solution of (6.226, c, /) for M is continuous across 
the boundaries of region IV. For vacuum solutions, it is appropriate to 
use (7.8) and equations (7.9) may then be written in the form 



S f = -2(f + g)(l-EE)- 2 E f E f 
S g = -2(f + g)(l-EE)- 2 E g E g . 



(11.21) 
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To ensure that the boundary conditions are satisfied, it is then essential 
that the solution of these equations should include the necessary compo- 
nents (7.10). 

In this approach, solutions in the interaction region for Z(t, z) or 
E(t, z) are related to the functions V(f, g) and W{f, g), and hence S may 
be obtained as a function of / and g. The boundary conditions described 
in Chapter 7 may then be considered to place restrictions on the structure 
of these functions, which characterize the approaching waves. 

In practice, these boundary conditions are difficult to apply, basically 
because the condition that M be continuous is only indirectly applied to 
the functions V and W, or Z, or E, that feature in the main equations. 
It is therefore often convenient to apply the boundary condition in the 
form (7.15) or (7.16). Writing Z = Z(u,v), this becomes 



lim 
lim 



ZyZ V 



n 2 (n 2 - l)v n *- 2 (Z + Z) 2 
2 Z U Z U 

m(m - i)u n i- 2 (Z + Z) 2 



cl 



(11.22) 



This form is particularly convenient when n\ = n 2 = 2, which occurs when 
impulsive waves are present. The boundary conditions in this limited case 
only, have been discussed by Ernst, Garcia-Diaz and Hauser (19876). 1 
In terms of the functions / and g, the conditions (11.22) become 



lim 

9-1/2 
/-1/2 

lim 

/-1/2 
9-1/2 



a -9) 



Z g Z g 


k 2 


(Z + Z) 2 _ 


2 


z f z f 


_ fc i 


(Z + Z) 2 _ 


2 



(11.23) 



where k\ and k 2 must satisfy the inequalities (7.13). 

Alternatively, writing E = E(u,v), the boundary conditions require 

that 

2 



hm 
lim 



EyEy 



n 2 (n 2 - l)u»3- 2 (1 - EE) 2 _ 
2 



E U E U 



ni(ni - l)w n i- 2 (1 - EE) 2 



= Cn 



cl 



(11.24) 



1 For an alternative approach to the formulation of the boundary conditions 
see Hauser and Ernst (1989a). 
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which, in terms of the functions / and g become 



lim 

/-1/2 

lim 

/-1/2 
9—1/2 



(§-<?); " " 



a-/) 



(1 — EE) 2 

E f E f_ 
(l-EE)< 



2 
2 



(11.25) 



Finally, it is convenient to express the boundary conditions in terms 
of the variables t and z. These become, for Z = Z(t, z), 




(11.26) 



and for E = E(t, z), 



lim 

2^0 

lim 

2^0 



(E t + E z )(E t + E z ) 
(1 - EE) 2 

(E t -E z )(E t -E z ) 
(I -EE) 2 



= 2k! 
= 2k 2 



(11.27) 



where k\ and k 2 satisfy (7.13). 



11.3 Colinear solutions 

In the next chapter, the approach described above will be used to derive 
new solutions for colliding gravitational waves whose polarization vectors 
are not aligned. Such solutions essentially have W non-zero, and hence Z 
and E are complex. However, before moving on to consider such cases, it 
is appropriate first to review the colinear solutions described in previous 
chapters. These solutions have Z and E real. 

(i) The solution of Khan and Penrose (1971), discussed in Chap- 
ter 3 and Section 8.2, which describes the collision of aligned impulsive 
gravitational waves, is given by 



(11.28) 
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In this case / = \ — u 2 and g = \ — v 2 . It may be observed that, for 
stationary axisymmetric space-times, (11.28) is the Ernst potential which 
leads to the Schwarzschild solution. 

(ii) The Szekeres (1972) class of solutions, described in Chapter 9, is 
given by 

/ 1 ,+\ (fci+feO/2 /-, , x (fci-fc 2 )/2 



1-tJ 

In this case / = \ — (au) ni + . . . and 9 = \ — (bv) n2 + . . ., with k 2 = 
2(1 - 1/m) for i = l,2 and m > 2. 

(iii) The 'solution' of Stoyanov (1979) given by (10.24) uses 

Z= (l-t 2 ) a/2 (l-z 2 ) a/2 . (11.30) 

In Section 10.2 it has been argued that this solution must be considered 
to be unphysical on its own, but it may be included as a factor in more 
general solutions. It does not satisfy the boundary conditions (11.26). 

(iv) The solution of Ferrari and Ibanez (1987a) and Griffiths (1987), 
which is described in Section 10.3, is characterized by 



1 + t 
1-t 



( n + t\ \ (1L31) 

where a = (3z 2 - 1) {^(3t 2 - 1) log \^±±) - |tj . 

In this case / = \ — {c\u) n + . . . and g = \ — (c 2 v ) n + . . ., where n > 2 
and (a + b) 2 = 2(1 — 1/n). It may be observed that, when a = 1, (11.31) is 
the Ernst potential which, for stationary axisymmetric solutions, leads to 
the solution of Erez and Rosen (1959) which describes the external field 
of a non-rotating body with a quadrupole moment. 

(v) The generalized solution of Ferrari and Ibanez (19876), described 
in Section 10.4, is characterized by 

Z = (j^j k (1 " t 2 ) a/2 (l - z 2 T 12 - (H-32) 

This can be seen to include the Stoyanov factor (11.30), and to reduce to a 
Szekeres solution with k\ = hi = k when a = 0. Again / = \ — (ciu) n + . . . 
and g = \ — (c2v) n + . . ., where n > 2 and k 2 = 2(1 — 1/n). The degenerate 
cases occur when k = 1 and a = ±1. 
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(vi) The solution of Tsoubelis and Wang (1989) given by (10.44) has 

Z = ( £±| J f J (1 - t 2)a/2 (1 _ (1L33) 

This can be seen to be a generalization of (11.29) and (11.32), and its 
properties can immediately be deduced. 

(vii) The odd order solution of Griffiths (1987) described in Section 
10.6 uses the Ernst potential 

Z = e - azQ ^\ where Q t (t) = l - log (jZ^j ~ 1 - ( n -34) 

Again / = \ — {c\u) n + . . . and g = \ — (c2v) n + . . ., where n > 2 and 
a 2 = 8(1 - 1/n). 

It is not difficult to see how further solutions of this type can be 
generated. 



